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CN ! Abstract 

The sum-capacity is studied for a iiT-user degraded Gaussian multiaccess relay channel (MARC) 
C/3 ' where the multiaccess signal received at the destination from the K sources and relay is a degraded 

version of the signal received at the relay from all sources, given the transmit signal at the relay. An outer 
bound on the capacity region is developed using cutset bounds. An achievable rate region is obtained 
\f^ ' for the decode-and-forward (DF) strategy. It is shown that for every choice of input distribution, the 
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rate regions for the inner (DF) and outer bounds are given by the intersection of two i^-dimensional 
polymatroids, one resulting from the multiaccess link at the relay and the other from that at the 
(<— ^ ■ destination. Although the inner and outer bound rate regions are not identical in general, for both 
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cases, a classical result on the intersection of two polymatroids is used to show that the intersection 
belongs to either the set of active cases or inactive cases, where the two bounds on the iiT-user sum- 

•i— I . 

^^ . rate are active or inactive, respectively. It is shown that DF achieves the capacity region for a class of 

jrt ■ degraded Gaussian MARCs in which the relay has a high SNR link to the destination relative to the 

multiaccess link from the sources to the relay. Otherwise, DF is shown to achieve the sum-capacity for 
an active class of degraded Gaussian MARCs for which the DF sum-rate is maximized by a polymatroid 
intersection belonging to the set of active cases. This class is shown to include the class of symmetric 
Gaussian MARCs where all users transmit at the same power. 
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Index Terms 

Multiple-access relay channel (MARC), outer bounds, achievable strategies, Gaussian and degraded 
Gaussian MARC. 

I. Introduction 

The multiaccess relay channel (MARC) is a network in which several users (sources) commu- 
nicate with a single destination in the presence of a relay [1]. The coding strategies developed for 
the relay channel [2], [3] extend readily to the MARC [4], [5]. For example, the strategy of [3, 
theorem 1], now often called decode-and-forward (DF), has a relay that decodes user messages 
before forwarding them to the destination [4], [5]. Similarly, the strategy in [3, theorem 6], now 
often called compress-and-forward (CF), has the relay quantize its output symbols and transmit 
the resulting quantized bits to the destination [5]. 

Capacity results for relay channels are known only for a few special cases such as the class 
of degraded relay channels [3] and its multi-relay generalization [6], [7], the class of semi- 
deterministic relay channels [8], the class of orthogonal relay channels [9], [10], the class of 
Gaussian relay without delay channels [11], [12], and the class of ergodic phase-fading relay 
channels [4]. For the class of degraded relay channels, the degradedness condition requires that 
the received signal at the destination be independent of the source signal when conditioned 
on the transmit and receive signals at the relay. For the Gaussian case, this simplifies to the 
requirement that the signal received at the destination be a noisier version of that received at 
the relay conditioned on the transmitted signal at the relay. This condition immediately suggests 
that requiring the relay to decode the source signals should be optimal. In fact, for this class, 
applying this degradedness condition simplifies the cut-set outer bounds to coincide with the DF 
bounds. For the MARC, we generalize this degradedness condition to requiring that the signal 
received at the destination be independent of all source signals conditioned on the transmit and 
receive signals at the relay. Applying this degradedness condition to the cutset outer bounds for 
a MARC, however, does not simplify the bounds to those achieved by DF. 

A K-u^QX Gaussian MARC is degraded when the multiaccess signal received at the destination 
from the K sources and relay is a noisier version of the signal received at the relay from all 
sources, given the transmit signal at the relay. For a K-u^qx degraded Gaussian MARC, we 
develop the DF rate region as an inner bound on the capacity region using Gaussian signaling 



at the sources and relay. The outer bounds on the capacity region are obtained by specializing 
the cut-set bounds of [13, Th. 14.10.1] to the case of independent sources [14] and by applying 
the degradedness condition. In fact, for each choice of input distribution, both the DF and the 
cutset rate regions are intersections of two multiaccess rate regions, one with the relay as the 
receiver and the other with the destination as the receiver. In general, however, the inner and 
outer bounds differ in their input distributions as well as the rate bounds. The outer bounds 
allow a more general dependence between the source and relay signals relative to DF where 
we use auxiliary random variables, one for each source, to relate the transmitted signals at the 
sources and relay. For the Gaussian degraded MARC, we show that using Gaussian input signals 
at the sources and relay maximizes the outer bounds. For the inner bounds, we use Gaussian 
signaling at the sources and the relay via K Gaussian auxiliary random variables. As a result, 
for each choice of the appropriate Gaussian input distribution, both the DF and outer bounds are 
then parametrized by K source-relay cross-correlation coefficients, i.e., a i^-length correlation 
vector. Specifically, each DF coefficient is a product of the two power fractions allocated for 
cooperation at the corresponding source and the relay, respectively. We show that the DF rate 
region over all feasible correlation vectors is a convex region. On the other hand, for the outer 
bounds, all the rate bounds at the relay except for the bound on the A'-user sum-rate are non- 
concave functions of the correlation coefficients, and thus, the outer bound rate region requires 
time-sharing. Finally, we also show that for every feasible choice of the correlation vector, the 
multiaccess regions achieved by the inner and outer bounds at each receiver are polymatroids, 
and the resulting region is an intersection of two polymatroids. 

We use a well-known result on the intersection of two polymatroids [15, chap. 46] to broadly 
classify polymatroid intersections into two categories, namely, the set of active and the set of 
inactive cases, depending on whether the constraints on the K-user sum-rate at both receivers 
are active or inactive, respectively. In fact, we use [15, chap. 46] to show that the K-user sum- 
rate for the inactive cases is always bounded by the minimum of the (inactive) K-user sum-rate 
bounds at each receiver, and thus, by the largest such bound. For both the inner and outer bounds, 
the intersection of the two rate polymatroids results in either an active or a inactive case for 
every choice of correlation vectors. In fact, the minimum of the K-user sum-rate bounds at the 
relay and destination is the effective sum-rate only if the polymatroid intersection is an active 
case and is strictly an upper bound for an inactive case. 



Irrespective of the above mentioned distinction, we first consider the problem of maximizing 
the minimum of the i^-user sum-rate bounds at the relay and destination over the set of all 
correlation coefficients. We solve this max-min optimization problem using techniques analogous 
to the classical minimax problem of detection theory [16, II.C]. We refer to a sum-rate optimal 
correlation vector as a max-min rule. 

For both the inner and outer bounds, we show that the max-min optimization described above 
has two unique solutions. The first solution is given by the maximum K-u^qx sum-rate achievable 
at the relay and results when the multiaccess link between the sources and the relay is the 
bottle-neck link. For this case, we show that the intersection of the rate regions at the relay and 
destination belongs to the set of active cases and is in fact the same as the region achieved at the 
relay. We further show that this region is the same for both the inner and outer bounds and is 
the capacity region for a class of degraded Gaussian MARCs where the source and relay powers 
satisfy the bottle-neck condition for this case. 

The second solution pertains to the case in which the bottle-neck condition described above 
is not satisfied, i.e., the fcT-user sum-rate at the relay is at least as large as that at the destination. 
For this case, we show that for both the inner and outer bounds the max-min optimization 
solution requires the i^-user sum-rate bounds at the relay and destination to be equal. In fact, 
we show that both the inner and outer bounds achieve the same maximum sum-rate for this case. 
Further, for both sets of bounds, we show that this maximum is achieved by a set of correlation 
vectors, i.e., the max-min rule is a set rather than a singleton. Recall, however, that the sum-rate 
computed thus is achievable for either bound only if there exists at least one max-min rule for 
which the polymatroid intersection belongs to the set of active cases; otherwise, the computed 
maximum is strictly an upper bound on the maximum sum-rate. Combining this with the fact that 
the maximum inner and outer fc'-user sum rate bounds for this case are the same, we establish 
that DF achieves the sum-capacity of an active class of degraded Gaussian MARCs, i.e., a class 
for which the maximum sum-rate is achieved because there exists at least one max-min rule for 
which the polymatroid intersection is an active case. We also show that the class of symmetric 
Gaussian MARCs, in which all sources have the same power, belongs to this active class. Finally, 
for the remaining inactive class of degraded Gaussian MARCs in which no active case results 
for any choice of the max-min rule, we provide a common upper bound on both the DF and the 
cutset sum-rates. 



This paper is organized as follows. In Section |ll] we present a model for a degraded Gaussian 
MARC. In Section UlI] we develop the cut- set bounds on the capacity region of a MARC. In 
Section UV] we determine the maximum K-user DF sum-rate. We discuss our results and conclude 
in Section |Vl 

II. Channel Model and Preliminaries 

A K-user degraded Gaussian MARC has K user (source) nodes, one relay node, and one 
destination node (see Fig. [T]). The sources emit the messages Wk, k = 1,2, . . . , K, which are 
statistically independent and take on values uniformly in the sets {1, 2, ... , M^}. The channel 
is used n times so that the rate of Wk is Rk = Bk/ n bits per channel use where Bk = log2 Mk 
bits. In each use of the channel, the input to the channel from source k is X^ while the relay's 
input is X,.. The channel outputs Yj. and Y^, respectively, at the relay and the destination are 

>;=(X]xj+z, (1) 

ld= (X^Xfcj +X, + Zrf (2) 

= Yr + Xr + Za (3) 

where Zr and Z^ are independent Gaussian random variables with zero means and variances 
Nr and N^, respectively, such that the noise variance at the destination is 

Nd = Nr + NA. (4) 

We assume that the relay operates in a full-duplex manner, i.e., it can transmit and receive 
simultaneously in the same bandwidth. Further, its input X.^ in each channel use is a causal 
function of its outputs from previous channel uses. We write /C = {1, 2, . . . , K} for the set of 
sources, T = /C U {r} for the set of transmitters, TZ = {r, d} for the set of receivers, Xg = {Xk 
: k E S} for all S C JC, and S'' to denote the complement of S in /C. 

The transmitted signals from source k and the relay have a per symbol power constraint 

E[\Xkf]<Pk keT. (5) 




Fig. 1. A two-user Gaussian degraded MARC. 



One can equivalently express the relationship between the input and output signals in ([3]) as a 
Markov chain 

(Xi,X2,...,X^)-(i;,X,)-F,. (6) 

For K = 1, ^ simplifies to the degradedness condition in [3, (10)] for the classic (single source) 
relay channel. A degraded Gaussian MARC is symmetric if Pk = P, for all k. Thus, a class of 
symmetric DG-MARCs is characterized by four parameters, namely, P, Pr , Nr, and A^^. 

The capacity region Cmarc is the closure of the set of rate tuples (-Ri, R2, . . . , Rk) for which 
the destination can, for sufficiently large n, decode the K source messages with an arbitrarily 
small positive error probability. As further notation, we write Rs = Xlfee^-^'^ ^^^ ^ ~ O^r, ^d)- 
We write and 1 to denote vectors whose entries are all zero and one, respectively, and C{x) = 
log(l + x)/2 to denote the capacity of an AWGN channel with signal-to-noise ratio (SNR) x. 
We use the usual notation for entropy and mutual information [13], [17] and take all logarithms 
to the base 2 so that in each channel use our rate units are bits. 



III. Outer Bounds 

An outer bound on the capacity region of a MARC is presented in [14] using the cut-set 
bounds in [13, Th. 14.10.1] as applied to the case of independent sources. We summarize the 
bounds below. 

Proposition 1: The capacity region Cmarc is contained in the union of the set of rate tuples 



(i?i, R2, . . . , Rr) that satisfy, for all 5 C /C, 



Rs<mm{I{Xs;Yr,Yd\Xs^,Xr,U),I{Xs,Xr;Ya\Xsc,U)} 



(7) 



where the union is over all distributions that factor as 



P{u) ■ \\\^^^P{Xk\u)j ■p{Xr\xKi,u) ■p{yr,yd\xK,Xr)- (8) 

Remark 1: The time-sharing random variable U ensures that the region in (|7]) is convex. 
One can apply Caratheodory's theorem [18] to this K-dimensional convex region to bound the 
cardinality of f/ as |W| < i^ + 1. 

Consider the outer bounds in Proposition [TJ For a degraded Gaussian MARC applying the 
degradness definition in ^ simplifies (|7]) as 



Rs < min {/(X5; Yr\XrXscU), /(X5X,; Ya\XsM)} for all 5 C /C 



(9) 



for the same joint distribution in ([8]). In the following theorem, we develop the bounds in Q 
with f/ as a constant. For notational convenience, for a constant U, we write fi^.^ and Bd^s to 
denote the first and second terms, respectively, of the minimum on the right-side of dH). The 
proof of the following theorem is detailed in Appendix HI 



Theorem 1: For a degraded Gaussian MARC, the bounds Br,s and 5^,5 are given by 
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where 7^, = 1 - Y.k&s- 7fc and 



(10) 



(11) 



^/^^iJ^P;. = E{XkXr) forallA;G/C. (12) 

Remark 2: For K = 1, the bounds in (flOl) and ([TT]) simplify to the first and second bound, 



respectively, for the degraded relay channel in [3, theorem 5]. 

Remark 3: The source-relay cross-correlation variables 7^, for all k, satisfy (11051) . i.e., they 
lie in the closed convex region Tqb given by 

OB= Sli^-'^lk<l> ■ (13) 

k&K. ) 

The bound 5^,5 in (flOl) . in general, is not a concave function of 7 for any S C JC. For a fixed 
7 , in Appendix HVl we show that fij.5 is a concave function of 7 . This in turn implies that 
Br,jc is a concave function of 7 . Further, in Appendix Hill we show that for all S, B^^s in (fTTI) 
is a concave function of 7 . 

Remark 4: In the expression for B^^s in (HTI) . the terms involving the cross-correlation coef- 
ficients quantify the coherent combining gains that result from choosing correlated source and 
relay signals. On the other hand, the expression for Br^s in (flOl) quantifies the upper bounds 
on the rate achievable at the relay when one or more source signals are correlated with the 
transmitted signal at the relay. 

The rate region TZqb enclosed by the cut-set outer bounds is obtained as follows. From (11191) 
for any choice of 7 , the rate region is an intersection of the regions enclosed by the bounds 
Br^s and Bd^s for all S. Since Br^s is not a concave function of 7 , one must also consider 
all possible convex combinations of 7 to obtain TZqb- For the fi'-dimensional convex region 
TZoB, one can apply Caratheodory's theorem [18] to express every rate tuple (-Ri, R2, ■ ■ ■ , Rr) 
in TZoB as a convex combination of at most K + 1 rate tuples, where each rate tuple is obtained 
for a specific choice of 7 . Let 9 denote the collection of all vectors 77 that satisfy 

Y:ltXVm = 1 (14) 

and let ( = ({7^}ii:+i,^) G Fq^^ x 6 denote a collection of i^ + 1 power fractions and weights 
such that the rate tuple achieved by the m*'^ vector 7^^ is weighted by the m*^ non-negative 
entry of the weight vector i], for all m E /CU {K +1}. Finally, since Tqb in (fT3l) is a closed 
convex set, J2m=iVmJ^^ £ ^ob- The following theorem presents an outer bound on the capacity 
region of the degraded Gaussian MARC. 

Theorem 2: The capacity region Cmarc of a degraded Gaussian MARC is contained in the 
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Fig. 2. Five possible intersections of TZr and TZd for a two-user Gaussian MARC. 



Case 3c (Active) R2 



region TZqb given as 



^oB= U (7^f (c)n7^^^(c)) 

Ceros 



where the rate region TZ'f (C) , j = r, d, is 



nf (c) = {{Ru i?2, . . . , Rk) ■.o<Rs< Bj,s (C) } 



and the bound Bj^s is given by 



K+l 



jS:' 



(15) 



(16) 



(17) 



B,,s (C) = E ^™^^-.^ (: 

m=l 

Theorem 3: The regions 7?.°'' (C) and 7^|^^ (C) are polymatroids. 

Proof: In Appendix HI] we show that for each choice of input distribution satisfying ([8]), the 
bounds in (ISTI) are submodular set functions, i.e., they enclose regions that are polymatroids. For 
the optimal Gaussian input distribution, this implies that 7^°^ ((^) and TZ^^ ((^) are polymatroids 
for every choice of ^. ■ 
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The region TZqb in (1611) is a union of the intersections of the regions 71° and 7^^ , where 
the union is taken over all convex combinations of 7 . Since TZqb is convex, we obtain the 
boundary of TZqb by maximizing the weighted sum ^k&f^^H-kRk over all Tqb and for all /ifc > 0. 
Specifically, we determine the sum-rate Rfc when fj.^ = 1 for all k. In general, to determine the 
intersecting polytope, one has to consider all possible polytope shapes for the regions TZ°^ and 
1Z°J'. However, since 71°^ and W^ are polymatroids, we use the following lemma on polymatroid 
intersections [15, p. 796, Cor. 46.1c] to broadly classify the intersection of two polymatroids into 
two categories. The first inactive set category includes all intersections for which the constraints 
on the two i^-user sum-rates are not active. This implies that no rate tuple on the sum-rate plane 
achieved at one of the receivers lies within or on the boundary of the rate region achieved at the 
other receiver. On the other hand, the intersections for which there exists at least one such rate 
tuple, i.e., the constraints on the two i^-user sum-rates are active in the final intersection, belong 
to the category of active set. In Fig. [21 for a two-user MARC we illustrate the five possible 
choices for the sum-rate resulting from an intersection of 7^°^(7^) and Tl^^('j ). Cases 1 and 
2 belong to the inactive set while cases 3a, 36, and 3c belong to the active set. We henceforth 
refer to members of the active and the inactive sets as active and inactive cases, respectively. 
Note that Fig. [2] illustrates two specific 7^°'' and Tl^f' polymatroids for cases 3a, 3b, and 3c. In 
general the active set includes all intersections that satisfy the definition for this set including 
cases such as TZ"'' C TZ^ and vice-versa. Finally, note that the sum-rate is a minimum of the 
sum-rates at the two receivers only for the active cases 3a, 36, and 3c. For the inactive cases 
1 and 2, the Ri + R2 constraints are no longer active and the sum-rate is given by the bounds 
Br,{2} + Bd,{i} and -Br,{2} + Bd,{i}, respectively. We use the following lemma on polymatroid 
intersections to generalize this observation and develop an outer bound on the K-user sum-rate. 

Lemma 1: Let Rs < /i (S) and Rs < /2 {<S), for all 5 C /C, be two polymatroids such that 
/i and /2 are nondecreasing submodular set functions on /C with /i (0) = /2 (0) = 0. Then 

m^xRic = min {fi (S) + /s (/C\5)) . (18) 

SQJC 

From Lemma[i|we see that the maximum i^-user sum-rate Rjc that results from the intersection 
of two polymatroids, Rs < /i (5) and R^ < /2 (S) is given by the minimum of the two K-user 
sum-rate planes /i (/C) and /2 (/C) only if both the sum-rates are at most as large as the sum of the 
orthogonal rate planes /i (S) and /2 (/C\5), for all 7^ 5 C /C. Further, the resulting intersection 
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belongs to the set of active cases. Conversely, when there exists at least one 7^ 5 C /C for 
which the above condition is not true, an inactive case results. Physically, an inactive case results 
when a subset S of all users achieve better rates at one of the receivers while the remaining 
subset of users achieve a better rate at the other receiver. For such inactive cases, the maximum 
sum-rate in (fTSi) is the sum of two orthogonal rate planes achieved by the two complementary 
subsets of users. As a result, the K-user sum-rate bounds /i(/C) and /2(/C) are no longer active 
for this case, and thus, the region of intersection is no longer a polymatroid with 2^^ — 1 faces. 

In the following theorem we use Lemma \T\ to develop the upper bound on the K-user sum- 
rate. For a Gaussian input distribution, the polymatroids TZ°^ and 7?.J^^ are parametrized by (^, 
and thus. Lemma [T] applies for each choice of (^. 

Theorem 4: For each ( G Tqb, the maximum K-user sum-rate Ric resulting from the inter- 
secting polymatroids T^"** and 7^^^ is 

{BdA + Br AC condition 1 

mm {Br,ic, Bd,K) otherwise 

where 5^,5 and Brs for all S are functions of ( and condition 1 is defined for any $ ^ A C JC 
as 

Bd,A + Br, A- < min {Br,K, Bd,jc) ■ (20) 

Remark 5: The condition in (|20|) determines whether the intersection of two polymatroids 
belongs to either the set of active or the set of inactive cases with respect to the A^-user sum- 
rate. 

Proof: The proof follows from applying Lemma [T] to the maximization of Rjc for each 
choice of C- ■ 

For a fixed transmit power P^, for all k E T, and noise variances A^,. and A^^, the choice 
of C determines whether the intersection of TZ"'' (C) ^nd T^J^'' (C) belongs to the set of active or 
inactive cases. For each choice of (, from Theorem |4] an active case results only if for all 2^—1 
non-empty subsets A of /C, the condition in (|20|) does not hold. Further, for any (^ that results 
in an inactive case, from Theorem HI the sum-rate is bounded as 

Bd,A + Br,A- < min {Br^jc, Bd,K) < max min {Br,jc, Bd,K.) ■ (21) 

Ceros 
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To this end, we consider the optimization problem 

Rk: = max min (5^,^ (() , Bd,ic (C)) • (22) 

In general, optimizing non-convex functions is not straightforward. However, since Bj-^/c and 
B(i^!C are concave functions of 7 , the above max-min optimization simplifies to 

Note that the optimization is performed over the same set in (l22l) and (|23]) as Fob is a closed 
convex set. In Appendix |Vl we show that the max-min problem in (l23l) is a dual of the classical 
minimax problem of detection theory, (see for e.g., [16, II.C]). This allows us to apply the 
techniques used to obtain a minimax solution to maximize the bounds in (|23l) over all 7 in 
Tqb (see also [9]). We write 7* to denote a sum-rate optimal allocation, i.e., a max-min rule, 
and write Q to denote the set of all 7* maximizing (|23] ). A general solution to the max-min 



optimization in (1231) simplifies to three cases [16, II.C]. The first two correspond to those in 
which the maximum achieved by one of the two functions is smaller than the other, while the 
third corresponds to the case in which the maximum results when the two functions are equal 
(see Fig. |4l). For Br^K and Bd,K defined in (flOl) and ([TT]) . respectively, we now show that the 
solution simplifies to the consideration of only two cases. The following theorem summarizes 
the solution to the max-min problem in (l23l) . The proof is developed in Appendix |Vl 
Theorem 5: The max-min optimization 

^^ ^ 7,'^roB "'''' {^'''' fc'c) ' B'^''^ [if^ ] (24) 

simplifies to the following two cases. 

Nr I ' 



Case 1: R^ = C \ ^^ I , B^,^ (0) < B^,^ 



(25) 



fee/c 



Case 2: i?^ = C E f - ^^^ ^ i?*, i?.V fc) = i?^,^ (l* 



where Pmax = maxfcgK:-Pfc, A^ = Pfe/Pmax, and x* = Y^keK V^kll is the unique solution 
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satisfying Brx f 7^ = Brx f 7^ and is given by 






with 



^ p V p (27) 

^2 - TV^ + IVd' ^"^ ^3 - W, • 

Remark 6: The maximization in (|24l) is independent of whether the optimal 7* results in an 
active or an inactive case. However, not all max-min rules 7* G ^ will result in an active case. 
In general, active cases may be achieved only by a subset Qa '!^ G- However, irrespective of 
the kind of intersection, from Lemma [H (l25l) is an upper bound on the i^-user sum-rate cutset 
bounds. 

In the following theorem we show that it suffices to consider two conditions in determining 
the largest outer bound on the i^-user sum-capacity. We enumerate the two conditions as 

Condition 1: Br,ic{0) < Bd,K{0) 

(28) 
Condition 2: 5,,,yc(0) > Ba,ic{0). 

The first condition implies that the maximum i^-user cutset bound at the relay is smaller than 
the corresponding bound at the destination; for this case, we show that B^^siQ) < Bd,siQ) f^^ 
all (S C /C, i.e., TZqb = T^"'' *^ T^d'- ^^ the other hand, when condition 2 occurs, i.e., when 
condition 1 does not hold in (|28l) . we use the monotone properties of Br^/c and Bd^/c and Lemma 
[U to show that 

with equality achieved in (|29l ) when the polymatroid intersection is an active case. From Theorem 
Owe have that a continuous set, G, of 7* maximizes the right-hand- side of (|29l ). We show that 



the bound in (1291 ) is achieved with equality when there exists a 7* that results in an active case, 
i.e., in a non-empty Qa- Finally, for the class of symmetric degraded G-MARCs, we prove the 
existence of an active case that maximizes the sum-rate. 



14 



Theorem 6: The largest outer bound W^ on the K-user sum-rate is 
Rt < C ((E^ft) - ^^T^) ' otherwise 



(30) 



where x* = J2keK V^^^'^l ^^ ^^^ unique solution satisfying Br^i^j*) = -Brf,/c(7*) and is given 
by (|26l) and (l27l) . The bound in (|30l) is achieved with equality only when the intersection of 



Tlf'{'y*) aiid Tlf{Y,c) results in an active case. The bound is achieved with equality for the class 
of symmetric degraded G-MARCs. 

Proof: Let 7^ = 0. From ^ we see that Bs ril ^ 7^ 0) < Bs r(0), for all 5 C /C, i.e, the 
region TZr (7r) i^ largest at 7 = 0. Expanding Bs,r and Bs^d at 7 = from (fTOl) and ([TT]) . 
respectively, we have 

B.,^ (0) = C (^^) (31) 

B.. (0) ^ C (^ + I) . (32) 

The sum-rate resulting from the intersection of 7^°^ (0) and TZ'd (0) falls into one of following 
two cases. 

Case 1: The first case results when Bjc-r (0) < Bfc,d (0). From (|3T1) and (l32l) this condition 
simplifies to 

Expanding (l33l) . we have, for any 5 C /C, 

EPk + Pr 

< '-^^ (35) 

where (l35l) follows from the degradedness condition in dH). Thus, -Br,K:(0) < -Bd,/c(0) implies that 
B,.s{0) < Bd,siO) for all S C JC, i.e., 7^°''(0) C 7^2^(0), and 7^oB(0) = 7^^^(0). The maximum 
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i^-user sum-rate upper bound for this active case is then 



Rt = Br,K = C{Y.keKPk/Nr). (36) 

Case 2: The second case results when B/c^r (0) > Bjc,d (0) , i.e., when 

^ Pk Yj Pk 

Unlike case 1, (1371 ) does not imply that B^^r (0) > B^^ (0) or vice-versa, for all 5 C /C. From 
Theorem |4l the intersection of 7?.°'' (0) and T?.^^ (0) can result in either an active or an inactive 
case and thus, from (l20l) . we have 

R'^ < min(5,,yc (0) , B.^ic (0)) = B,,jc (0) (38) 

with equality for the active case. Note that from symmetry an active case results for the symmetric 
G-MARC. We now show that the sum-rate is increased for a 7* 7^ such that Br^jc il*jc) ~ 
Bd,jc [l*fc)- To simplify the exposition, we write Br^jc and B^x in (flOl) and (fTTj) as 

Pr,^(^) = C\'-^-^^-^] (39) 

S.,.M=C|^^ + -+ ^__ I (40) 



where 



A 
X = 



5^ Vt^ (41) 



fe=i 



and Afc = Pfe/Pmax where P^ax = max^e/c ^fc, for all /c. For all 7^ G [0, 1], we have 



|2L = ^ A;G/C (42) 

|!f = -^ A; e /C (43) 

7T^ = fcy^J. (44) 
97fe97j ' -^ 
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Thus, X is a concave function of 7 over the hyper-cube 7^ G [0, 1], for all k, and therefore, 
is concave for all 7^ satisfying (fT3l) . Further, from (fT3l) . we see that x is maximized when the 
entries of 7 satisfy J2k=i'yk = 1- Using techniques similar to those in Appendix HIH one can 



show that X achieves its maximum for a 7' with entries 



7l 



TK — 



Efc=l ^k 



for all k. 



and thus, we have 



X G 



' V eL 



c [o,Vir]. 



(45) 



(46) 



The functions Br^K (a^) and -B^^/c (a^) in (l39l) and (l40l) are monotonically decreasing and increasing 
functions of x, respectively. Substituting (l45l) in (flOl) . we have i?r,5(7'^) = for ^U S C IC. 
Thus, for the case in which i?/c,r (0) > -B/c,d (0), one can shrink the region TZ°'' from TZf' (0) 
just sufficiently such that for some 7* 7^ 0, _B,.,/c (7r) ^ -^'^.'c (^r) > -^'c.d (0)- In Theorem [6] 
we show that B/c^r = Bjc^ is maximized by a set Q of 7* satisfying 



Ik 



E7^A, 



X 



A:G/C 



(47) 



where x* is the unique value satisfying the quadratic Bjc,r{x) = Bjc,d{x). For 7* 7^ 0, from 
CI) one can verify that Br,s (7;^) < 5^,5 (0) for all 5, i.e., TZf (y^ C 7^f (0). On the other 
hand, substituting 7* in (fTTI) . 5^ 5 for all 5 C /C simplifies to 



-L/C 



B. 



d,s 7 



^yc 



C 



keS k£S 



(48) 



Comparing S^ 5 (0) in (|32l) with i?^ .5 ( 7t ) in (l48l) above, one cannot in general show that 



7^ 



ob 



%C 



D TZ'^ (0). In fact, the 7* chosen will determine the relationship between Bd,s{lZ 



-LJC 



-J-K' 



ob 



and -8^,5(0) for any S. Thus, for any 7* that equalizes fir^/c and Bd,K the polytope 7^°'' fl 71 
belongs to either the set of active or inactive cases. Recall that we write Qa to denote the set 



of 7* that results in an active case, i.e., the set of 7* for which the condition in (|20l) does not 
hold for all 2^ — 1 non-empty subsets A of fC. From Theorem IH we have that the sum-rate for 
the inactive case is always bounded by the maximum sum-rate developed in Theorem |5l Thus, 
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the maximum K-user sum-rate when -Br,A:(0) > Bd,/c{0) is 

Rr = < — — (49) 

^ maxB,^AO + BrMO<B* Ga = ^ 

where B* is defined in Theorem |51 We now show that for the class of symmetric G-MARC 
channels the bound B* is achieved, i.e., Qa 7^ 0- For this class since Pk = P for all k E }C, from 
symmetry B* can be achieved by choosing 7^ = 7* for all k such that from (|4TI) . we have 



Y = {x*f/K\ (50) 

From (|46|) . since < x* < ^/K, there exists an 7* < 1. From (fT3l) . we also require 7* < 1/K. In 
Theorem [T2] in Section HVl below, we prove the existence of a 7* < l/K for symmetric channels. 
From symmetry, since no subset of users can achieve better rates at one receiver than the other, 
the resulting 7?.,. (7*) fl TZd (7*) belongs to the set of inactive cases. The K-user sum-rate cutset 
bound for this class is given by the B* in (l25l) with Pmax = P and A^ = 1 for all k E JC. 
Finally, from continuity, one can expect that for small perturbations of user powers from the 
symmetric case, an active case will result. However, for arbitrary user powers, it is possible that 
Qa = 05 i-e-5 the set of all feasible 7* results in non-inactive cases. In general, however, obtaining 
a closed-form expression for the maximum sum-rate for the inactive cases is not straightforward. 



IV. Decode-and-Forward 

A DF code construction for a discrete memoryless MARC using block Markov encoding and 
backward decoding is developed in [4, Appendix A] (see also [19]) and we extend it here to the 
degraded Gaussian MARC. We first summarize the rate region achieved by DF below. 

Proposition 2: The DF rate region is the union of the set of rate tuples {Ri, R2, ■ ■ ■ , Rk) that 
satisfy, for all S C }C, 

Rs < min {I{Xs; Yr\XscV,cXrU), /(X^X,; Yd\XscVsM)} (51) 

where the union is over all distributions that factor as 

P(") ■ (llk=iPi'"k\u)p{xk\vk,u)j -pixrlvicu) ■ p{yr,yd\xr) ■ (52) 
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Proof: See [19]. ■ 

Remark 7: The time-sharing random variable U ensures that the region of Theorem [21 is 
convex. 

Remark 8: The independent auxiliary random variables V^, A; = 1, 2, . . . , i^, help the sources 
cooperate with the relay. 

For the degraded Gaussian MARC, we employ the following code construction. We generate 
zero-mean, unit variance, independent and identically distributed (i.i.d.) Gaussian random vari- 
ables Vfc, Vfc 0, and Vr^, for all k E JC, such that the channel inputs from source k and the relay 
are 

Xk = V^J\Vk,o + y^il- ak)PkVk, keIC, (53) 



K / K 



Xr=E VW^Vk + W 1 - EPk PrVr,0 (54) 

fc=i y V k=i J 

where a^ G [0, 1] and I3k G [0, 1] are power fractions for all k. We write 

Q-fc = ("1, "2, •••, aK] (55) 



§_^={Pi. /32, ..., Pk) (56) 

and 

r = { («yc,^^) : «fc e [0, 1], < Efce^c A < 1 for all k. } (57) 

for the set of feasible power fractions a^ and /? . Substituting (|53l) and (|54l) in (|5T1) . for any 
(tt/C'/^r) £ r» we obtain 

i?5<min(^/,,5(«,c).^rf,5(a/c.^yc)) for all 5 C /C (58) 
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where Ir^s and I^s, the bounds at the relay and destination respectively, are 







2EA/(l-«fc)A 



k- 



Nd kes\ NdNd 



(59) 



(60) 



From the concavity of the log function it follows that J^.^, for all S, is a concave function of 
ay^. In Appendix Hill we show that Id^s is a concave function of a^^ and P . The DF rate region, 
TZdf, achieved over all {af^, (3 ) G F, is then given by the following theorem. 
Theorem 7: The DF rate region IZdf for a degraded Gaussian MARC is 

nDF= U (7^r(«J^7^, (a^,^^)) (6i) 

where the rate region TZt, t = r,d, is 

^t(«yc,^J = { {RuR2,....RK).^<Rs<ks[a^.^^, forall5C/C }. (62) 

Proof: The rate region TZdf follows directly from Proposition [2l the code construction in 

(t53l) -(l54l). and the fact that Ir,s and Id,s are concave functions of {a,^,P^). ■ 

Theorem 8: The rate region TZdf is convex. 

Proof: To show that TZdf is convex, it suffices to show that Ir^s and Id,s, for all S, are 

concave functions over the convex set F of (a/c/^ir)- This is because the concavity of Ir,s and 

Id^s, for all S, ensures that a convex sum of two or more rate tuples in TZdf, each corresponding 

to a different value of {a,^, (3 ) tuple, also belongs to TZdf, ^■^■, satisfies (|62l ) for t = r,d. ■ 

Theorem 9: The rate regions TZr and TZd are polymatroids. 

Proof: In Appendix |II] we show that for every choice of input distribution satisfying (|52l) 

the bounds in (|5TI) are submodular set functions, and thus, enclose regions that are polymatroids. 

For the Gaussian input distribution in (l53l) and (l54l) . this implies that 7^^. («) and TZd (a, /?) are 

polymatroids for every choice of (a, /3), i.e., TZr and 7?.^ are completely defined by the comer 

(vertex) points on their dominant i^-user sum-rate face [15, Chap. 44]. ■ 
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The region TZof in (|6Tl) is a union of the intersection of the regions TZr and TZd achieved 
at the relay and destination respectively, where the union is over all {a^,l3) E T. Since TZdf 
is convex, each point on the boundary of TZdf is obtained by maximizing the weighted sum 
^keK.^^kRk over all F, and for all /i^ > 0. Specifically, we determine the optimal policy {a^, /?* ) 
that maximizes the sum-rate Ric when fik = 1 for all k. From (|6T1) . we see that every point on 
the boundary of TZdf results from the intersection of the polymatroids TZr{a,(^) and TZd^a^, (3 ) 
for some {a^,l3). Since TZr and TZd are polymatroids, as with the outer bound analysis, here 
too we use Lemma [T] on polymatroid intersections to broadly classify the intersection of two 
polymatroids into the categories of active and inactive sets. In the following theorem we use 
Lemma [T] to write the bound on the K-u^qx DF sum-rate. We remark that TZr and TZd are 
polymatroids parametrized by {a^, f3 ), and thus. Lemma [U applies for each choice of (a^^, f3 ). 

Theorem 10: For any {ax:,/3), the maximum i^-user sum-rate Rjc resulting from the inter- 
secting polymatroids 7?., and TZd is 

{TdA + TrA'^, condition 2 

min [Irx, h,K) , otherwise 

where condition 2 is defined for a 7^ ^ C /C as 

Td,A + IrA" < niin {Ir,K, Td,ic) ■ (64) 

Remark 9: The condition in (|64|) determines whether the intersection of two polymatroids 
belongs to either the set of active or inactive cases with respect to the i^-user sum-rate. 

Proof: The proof follows from applying Lemma [T] to the maximization Rj^ = Tlik&K^k 
for each choice of (ay^,/5). ■ 

We seek to determine the maximum sum-rate Rfc over all (a^^, /? ) G F. To this end, we first 
consider the optimization problem 

Rk. = ^ max min ( J^,^ (a,^) , Id,ic (a^, P^)) ■ (65) 

We write (aj^, /3*) to denote the max-min rule optimizing (|65l) and write P to denote the set of all 
(aj^,/9* ) maximizing (l23l) . A general solution to the max-min optimization in (|23l) simplifies to 
three cases [16, ILC]. The first two correspond to those in which the maximum achieved by one 
of the two functions is smaller than the other, while the third corresponds to the case in which 
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the maximum results when the two functions are equal (see Fig. Hj). For Ir^jc and Id^jc defined in 
(|59| ) and (l60l) . respectively, we can show that the solution simplifies to the consideration of only 
two cases. The following theorem summarizes the solution to the max-min problem in (l65l) . The 
proof is developed in Appendix |Vl 

Theorem 11: The max-min optimization 

Ric = max min (lr,K ("/c) ' ^dx (aye f^y-)) (66) 

simplifies to the following two cases. 



Case 1: Rk = C \ ^^ I Ir,K (i) < Id,jc (1, 0) 



K, + y/Kf + iK,~K2)Ko 



with 



^2 - ^, + TvT, and K3 - — ^- — . 



The entries of the optimal /3* are given by 



:i-"I.Ffc 






a;^ = i 



(67) 



Case 2: i?^ = C (( E^f ) - (^:^) ^ J* /.^ = ^.V 
where I*,^ = It,K:{Q$:, fT^), t = r,d, P^ax = max^ Pk with Afc = Pfc/Pmax, and 

g*^5^A,(l-«^) (68) 

is the unique value satisfying the quadratic Ir,K.{^K., P*,^) ~ ^d,K.{^K.i P*,^) and is given by 



(69) 



(70) 



PI = ^ l.fe=ili-"fejn. foj. all kE]C. (71) 



Remark 10: The optimal q* in (|69l ) is the same as that for the optimal x* in (1261) . Thus, from 
(l25l) and (|67|) . we see that for both cases, the maximum cutset bound is equal to the maximum 
DF bound on Rfc. 

From Lemma \\\ we see that the maximum sum-rate can be achieved by either an active or an 
inactive case. In the following theorem we show that it suffices to consider two conditions in 



22 



determining the maximum K-user DF sum-rate. We enumerate the two conditions as 

Condition 1: /.,yc(l) < V(1,0) 
Condition 2: Ir,ic (1) > Id,ic (i, 0) . 

The first condition implies that the maximum sum-rate at the relay is smaller than the corre- 
sponding rate at the destination; for this case, we show that Ir,s{X) < ^d,s{l,Q) for all S C JC, 
i.e., TZdf = T^r C TZd- Physically, this corresponds to the case where the relay has a high SNR 
link to the destination and the multiaccess link from the sources to the relay is the bottleneck 
link. Under this condition, we show that the sum-capacity of a degraded Gaussian MARC is 
achieved by DF. On the other hand, when condition 2 occurs, i.e., when condition 1 does not 
hold in (|72|) . we use the monotone properties of Ir^jc and Idx and Lemma [T] to show that 

Rk < max min { J,,,^ (a,^) , J^,^ ( a^, (3 )\ (73) 

with equality when the intersection of Tlr^a^) and ^^^^{a^.P) results in an active case. From 
Theorem [m a continuous set V of (aj^, (3* ) with a unique (3* for each choice of aj^ maximizes 
the right-side of (1731) . Furthermore, we show that the bound in (fTSl) is the sum-capacity when 
an active case achieves the maximum sum-rate. Finally, for the class of symmetric degraded 
G-MARCs, we prove the existence of an active case that achieves the sum-capacity. 
Theorem 12: The /C-user DF sum-rate Rk. for a degraded Gaussian MARC is 

For /r./cd) < Id,K.il,Q.), DF achieves the capacity region and the sum-capacity of the degraded 
Gaussian MARC. The upper bound on R/c in (1741) is achieved with equality only for a class of 
active degraded Gaussian MARCs for which there exists a {a*^, P*) E V such that 7?.r(aj^) fl 
T^di^lK^ (^*ir) ^^ an active case and is the sum-capacity for this class. This active class also includes 
the class of symmetric degraded Gaussian MARCs. 

Proof: Let a^ = l and (3 = 0. From (|59| ) and (l60l) . we see that J^,. and Is4 are 
monotonically increasing and decreasing functions of a^, respectively, for a fixed (3 , i.e., for 
any a^-* and a^-* satisfying (|57|) . with entries a]^' < o\' for all k e IC, TZr{(ljc ) — '^r(«AC ) 
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„(1) 



„(2) 



and TZdio^jc , P ,(■) ^ 'Tldi^cr^ , [3 ) . Thus, TZriai^) achieves its largest region for a^^ = I. The 
bounds /r,5 and Id,s can be expanded for this case using (|59| ) and (l60l) . respectively, as 



/r.5 



-^d,5 = C 









(75) 
(76) 



The resulting sum-rate satisfies one of two conditions and we enumerate them below. 

Condition 1: The first condition is /r,/c(i) < Id,K:il,Q)- From (l75l) and (1761) . this case requires 



E A E ^fc 



fce/c 



AT, 



Expanding (1771) . we have, for any S C JC, 



< 



keK 



N. 



El 



(77) 



T.Pk EPk + Pr EPkiNd-Nr] 



fce5 



K 



< 



< 



kes 



keS'^ 



Nd 

EPk + Pr 

kes 



NdNr 



N, 



(78) 



where ^ follows from ©. Thus, Ir,jc{l) < hxiLQ) implies that Ir,s{l) < Id,s{hQ.) for all 
S C IC, i.e., 7^r(l) C 7^d(l) and thus, 7^DF(1) = 7^^(l). Further, since 7^r(l) n 7^d(l,0) = 
7^r(l), the polymatroid intersection for this condition belongs to the intersecting set. Finally, 
recall that we chose /5 = 0. From (|59| ). we see that the choice of P does not affect TZr- Further, 
a non-zero 8^ does not increase Idtc- However, it can decrease Ids for some or all iS C /C as 



// 



'd,5 



hP. 



c 



\keS J V fcg^'^ 



\ 



</d,5(l,0) 



(79) 



/ 



thereby potentially decreasing TI-df (!)• Thus, for the condition in (1771) and from Theorem [51 
the i^-user sum-capacity of a degraded G-MARC for this case is 



Rk = Ir,K (1) = Br,^ (0) = C{Y.,^^Pk /Nr). 



(80) 



The max-min rule for this condition is (a*^, P*) = (1,0)- Finally, from condition 1 in Theorem 
[6] for a class of degraded Gaussian MARCs where the source and relay powers satisfy (|771) . DF 
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achieves the capacity region since 

nDF = nr{i) = nf{<d). (8i) 

Condition 2: The second condition requires //c,r(i) > -^yc,d(i,0), i.e., 

Yj Pk Yj Pk 

keic — k&c — -f^_ .g2) 

Unlike condition 1, one cannot show here that J^,, > I^^ for all 5 C /C or vice-versa. Thus, 
from Theorem \Wi the intersection of TZr (i) and 7?.^ (1 , 0) can result in either an active or an 
inactive case. From (l63l) in Theorem \Wi we then have 

R,c < min {Ir,K (1) , Id,jc (1, 0) } = Id,jc ( i, 0) (83) 

with equality for the active case. Note that from symmetry an active case results for the symmetric 
G-MARC. However, the bound on the sum-rate, and thus, the sum-rate too, can be increased 
using the fact that Ir^jc and Id,K: are monotonically increasing and decreasing functions of a^^, 
respectively. In fact, from (l59l) and (l60l) . we see that reducing some or all of the entries of a,^ 
from their maximum value of 1 reduces /,. y^ and either reduces or keeps unchanged some or all 
Ir^s while increasing I^^jc. Further, since Ir,siQ.) = for all 5 C /C, one can shrink the region TZr 
just sufficiently to ensure that there exists some a^ and /3* such that /r,K:(«Jc) = -^rf,K:(«5c'/^t-)- 
From Theorem nn Lr = Id,jc is maximized by a set V of (a*^, P*) where a*^ and P* satisfy 
(1681 ) and (TTTI ), respectively. Evaluating /^^ at a max-min rule (aj^,/^* ), we have 

EPk 



k es Ekes (^ - ^l) PkPr {l~al) PkPr 



h,s=c\^-^+ ^^-'i: , :r' " +^E ' \% - " i • (84) 



For a*^ 7^ 1, since /r,5, for all S, is a monotonically decreasing function of a^^ we have 
TZr (ajc) '^ ^r (D- On the other hand, comparing (1761) and (f84l) one cannot in general show 
that 7?.d (aj^,/9* ) ^ '^^(1,0). In fact, the aj^ chosen will determine the relationship between 
h,s{o^K.^ f^,^ and Id,s{h^ for any 5. Thus, for any (a;^,/5j^) that equalizes /,./c and J^^/c, the 
polytope IZr n T^rf belongs to either the set of active or inactive cases. Let Va ^V denote the set 
of (ttJc'/^ir) '^hat result in active cases. From Theorem [TOl we can write the maximum ii'-user 
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DF sum-rate when Ir,K.iV) > ^d,K.ihQ) as 

^'^='^ max h,Ai(^hP*) + IrM^h) < I\ ^a = ^^^^ 

where /* is given by (|67]) in Theorem [TTJ Finally, as shown in remark [TOl /* = B* where B* 
is the maximum outer bound sum-rate. 

We now show that for class of symmetric G-MARC channels, when the condition in (|82l ) holds, 
we achieve the /('-user sum-capacity. For this class, since P^ = P, from symmetry, J^x: = Ir,K. 
in (l60l) can be maximized by choosing al = a* for all k in (l68l) such that 



{l-a*) = {q*f/K. (86) 

From (t68l) . since < (g*)^ < X]^i -^fc = -^' ^h^^e exists an < a* < 1 that achieves /* in 
(1851) . Further, from symmetry, no subset of users achieves a larger rate at one of the receiver 
than any other subset, i.e., for al = a* and f3k = l/K, for all k, TZr fl TZd belongs to the set 
of active cases and the maximum K-user sum-rate for this class is P = B* . Recall that for the 
outer bound in Theorem [6l we need to prove that 7* G Vqb where 7* has entries 7* given by 



(1501) for all k. From (fT3|) and (1571) . we can write 

7, = (1 - au) h where (a^, ^^) G V. (87) 

. We then have 

J]7fc = J](l-afc)/3fc<l (88) 

k&K. fce/C 

where ([88]) follows from ([57]) and the fact that (1 - Ofc) A < A for all {a^,§_i^ G T. For the 
symmetric case, this implies that there exists a 7* = (1 — a*) /K satisfying (f88l) . In fact, for a* 
in ((86]), we obtain 7* = {q'f /K'^ = (x*)'^ /K^ < 1, i.e., the symmetric 7* in (|50l) is feasible 
and results in an active case. Since an active case achieves the same maximum sum-rate for both 
the inner and outer bound, we see that DF achieves the sum-capacity for the class of symmetric 
Gaussian MARCs. 

For the general case of arbitrary P^, from (f85l) and (|49l ) we see that DF achieves the maximum 
i^-user sum-rate outer bounds for an active class of degraded Gaussian MARCs for which 
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Rr{(lK:)^Rdi^K:, Pt-) belongs to the set of active cases. Further, DF achieves the same maximum 
value for all {a'^,P*) E Va ^ ^- In Appendix IVIIl we show that for the same choice of the 
K source-relay correlation coefficients for both the inner and outer bounds, the outer cutset 
bounds are at least as large as the inner DF bounds for all 5 C /C. This implies that for every 
(aj^,^) G Va, there exists a 7* with entries 

Y, = il-al)/3l for all A; (89) 

that results in an active case for the outer bounds, i.e., DF achieves the sum-capacity for the 
active class. Note that the outer bounds may also be maximized by other (a^^,/? ) that do not 
maximize the K-user DF sum-rate. 

Finally, as with the outer bounds, the optimization in (l85l) for P^ = is not straightforward. 
Further, comparing the DF and cutset bounds in (l85l) and (|49l) , respectively for the inactive 
cases, we see that the expression for the outer bounds involves time-sharing and can in general 
be larger than the DF bound. ■ 

It is straightforward to find numerical examples for condition 1 in Theorem [12] where DF 
achieves the capacity region. We focus on condition 2 and present two examples where DF 
achieves the sum-capacity of a two-user degraded Gaussian MARCs, with Va = V for one and 
Va CV for the other. 

Example 1: Consider a two-user degraded Gaussian MARC with Pi/Nj. = 6, P2/Nr = 4, 
Pi/Nd = 3, P2/Nd = 2, and Pr/N^ = 2. These SNR values satisfy the condition 2 given by 
(|82|) in Theorem [12] and thus, the DF sum-rate is maximized by a set of (a*^, P*) where a*^ 
satisfies 

(1 - a*) + ^ (1 - a;) = {q*f = 0.408, (90) 

and for every choice of a^ satisfying ([90] ). /3* is given by ([TT]) . The set of feasible a*^ has entries 
al e (0.83, 1] with a^ for each such al satisfying ([90]) such that a^ G (0.75, 1]. For these SNR 
parameters, the set Va = V and for each (aj^, /?* ) eV, the correlation values 7^ = (1 — a^) /5^, 
for all /c = 1, 2. result in the vector 7* G Qa- 

Example 2: We next consider a two-user example with Pi/Nj. = 6, P2/Nr = 0.4, Pi/Nd = 3, 
P2/Nd = 0.2, and Pr/Na = 2. These SNR values also satisfy the condition 2 given by ([82]) in 
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Theorem [T2l and thus, the DF sum-rate is maximized by a set of («5^,/3* ) where aj^ satisfies 



'—K.' 



(1 - a*) + ^ (1 - «;) = (g*)' = 0.197. (91) 

The set of feasible aj^ has entries a* G (0.96, 1] with ^2 for ^^^h such a* satisfying (|9T| ) such 
that a;2 ^ (0.416, 1]. Note that subject to (|9T| ). 0:2 decreases as ai increases and vice- versa. 
For these SNR parameters, the set Va consists of (a*^, /5* ) where the entries a^ and ^2 ^^ 
restricted to the sets (0.961, 0.979] and (0.731, 1], respectively. The remaining values for a* and 
0:2 satisfying (|9T1) result in a polymatroid intersection that belongs to the set of inactive cases. 
In fact, all such values result in the inactive case 2 illustrated in Fig. [2] for K = 2. 

Finally, for the two-user degraded Gaussian MARC, a numerical example illustrating Va = ^ 
does not appear straightforward despite using a wide range of ratios of Pi to P2, i.e., not all 
rate-maximizing intersections are such that one of the sources achieve better rates at one of the 
receivers while the other source achieves a better rate at the other receiver. A possible reason 
for this is because, at any receiver, the noise seen by both sources is the same, and thus, the 
source with smaller power typically achieves smaller rates at both receivers. It may be possible 
to increase the rate achieved at the destination by increasing the relay power; however, large 
values of relay power will result in the bottle-neck case where condition 1 in Theorem [T2l holds . 
Thus, it appears that it may always be possible to find an active case, particularly, one that 
maximizes the sum-rate. If this is true for any arbitrary K, then DF achieves the sum-capacity 
of the degraded Gaussian MARC. 

Remark 11: In the above analysis, we determined the sum-capacity for a degraded Gaussian 
MARC under a per symbol transmit power constraint at the sources and relay. One can also 
consider an average power constraint at every transmitter. The achievable strategy remains 
unchanged; for the converse we start with the convex sums of the outer bounds in (|7]) over 
n channel uses. In the z*^ channel use, the bounds at the relay and destination are given by 5^,5 
and Bfis in (flOl) and (fTTI) . respectively, for all S, except now the correlation parameters and 
power parameters are indexed by i. Recall that B^.s is a concave function of the correlation 
coefficients and power. On the other hand, Br^s for all 5 C /C is not a concave function of the 
power and cross-correlation parameters. However, we can use the concavity of B^^jc to show 
that the maximum bounds on the sum-rate in Thereom|4] remain unchanged. This in conjunction 
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with the steps in Theorem [6l lead to the same sum-capacity results. Finally, we note that as with 
the symbol power constraint, here too we require time-sharing to develop the outer bound rate 
region. 

V. Concluding Remarks 

In this paper, we have studied the sum-capacity of degraded Gaussian MARCs. In particular, 
we have developed the rate regions for the achievable strategy of DF and the cutset outer bounds. 
The outer bounds have been obtained using cut-set bounds for the case of independent sources 
and have been shown to be maximized by Gaussian signaling at the sources and relay. 

We have also shown that, in general, the rate regions achieved by the inner and outer bounds 
are not the same. This difference is due to the fact that the input distributions and the rate 
expressions for the inner and outer bounds are not exactly the same. In fact, the input distribution 
for the inner bound uses auxiliary random variables to model the correlation between the inputs 
at the sources and the relay and is more restrictive than the distribution for the outer bound. 
Despite these differences, in both cases the input distributions can be quantified by a set of K 
source-relay cross-correlation coefficients. Further, in both cases, we have shown that the rate 
region for every choice of the appropriate input distribution is an intersection of polymatroids. 
We have used the properties of polymatroid intersections to show that for both the inner and 
outer bounds the largest K-user sum-rate is at most the maximum of the minimum of the two 
K-user sum rate bounds, with equality only when the polymatroid intersections belongs to the 
set of active cases in which the i^-user sum rate planes are active. 

For both DF and the outer bounds, we have shown that the largest fi'-user sum-rate can be 
determined using max-min optimization techniques. In fact, we have shown that for both the 
inner and outer bounds the max-min optimization problem results in one of two unique solutions. 
The first solution results when the largest sum-rate from the K sources to the relay is the bottle- 
neck rate and for this case, we have shown that DF achieves the capacity region. We have further 
shown that the sum-rate maximizing polymatroid intersection for this case belongs to the set 
of active cases. Specifically, the sum-capacity as well as the entire capacity region is achieved 
by a max-min rule where the sources and the relay do not allocate any power to cooperatively 
achieving coherent combining gains at the destination, i.e., the auxiliary random variables Vk = 0, 
for all k. Thus, under Gaussian signaling, the capacity region is achieved by DF because the 
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inner and outer bounds at the relay, for V/c = 0, are I(Xs]Yr.\XrXsc) = I(Xs]Yr.\Xj.Xs':Vjc) 
for all 5 C /C (see © and dSB)- 

The second solution results when the largest sum-rates at the relay and the destination are 
equal. For this case, we have shown that DF achieves the sum-capacity for a class of active 
degraded Gaussian MARCs in which the sum-rate maximizing polymatroid intersection belongs 
to the set of active cases. We have also shown that this class of active degraded Gaussian MARCs 
contains the class of symmetric Gaussian MARCs. In general, for this class, we have shown that 
the max-min DF rule is such that V^ ^ for all k, i.e., a non-empty subset of sources and the 
relay divide their transmit power to achieve cooperative combining gains at the destination. We 
have also shown that the largest DF sum-rate is achieved by a relay power policy that maximizes 
the cooperative gains achieved at the destination, i.e., X^ is a unique weighted sum of Vk for 
all k where the weight for each source k is proportional to the power allocated by source k 
to cooperating with the relay. Our analysis has also shown that the maximum sum-rate admits 
several solutions for the power fractions allocated at the sources for cooperation subject to a 
constraint that results from the equating the two bounds on the sum-rate. For the outer bounds, 
we have shown that the K-user sum-rate outer bound is maximized by a set of cross-correlation 
coefficients that are subject to the same constraint as DF and the maximum sum-rate is the same 
as that for DF. Furthermore, for the class of active degraded Gaussian MARCs, we have shown 
that the set of DF max-min rules (a*^, P*) also maximizes the outer bounds by using the fact 
that the inner and outer bound coefficients can be related as 7^ = (l — ak) Pk, for all k. Finally, 
since a DF max-min rule requires a unique correlation between Xr and Vyc, conditioning the 
outer bound that uses Yr on Xr alone suffices to obtain the sum-capacity. 
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Appendix I 
Outer Bounds: Proof 

We now develop the proof for Theorem \T\ Recall that we write Br^s and Bd,s to denote, 
respectively, the first and second bound on Rg in @ for a constant U. Expanding the bounds 
on Rs in ^ for a constant U, we have 



Rs < mm {h{Yr\XrXs.) - h{Zr), h{Yd\Xs.) - h{Z, 



)}• 



(92) 



For a fixed covariance matrix of the input random variables X/c and Xr, one can apply a condi 
tional entropy maximization theorem [20, Lemma 1] to show that h{Yr\XrXsc) and h{Yd\Xs':] 
are maximized by choosing the distribution in ([8]) as jointly Gaussian. Consider the bound B^^s 
Expanding Yr, we have 

' E [var {Zkes^klXrXsc)] 



Rs<C 



Nr 



For Gaussian signals, using the chain rule, we have 



E 



var { Y^XklXrXsc 
kes 



detjKAic) 
det{KB\c) 



where 



(93) 



(94) 



A = 

B = [Xr] 

C = [Xs 



Y.k<^S^k Xr 



(95) 
(96) 
(97) 



and for random vectors X_ and Y_, the conditional covariance Kx\y_ is 



K 



K\Y_ 



E 



{X-E[X\Y]){X-E[X\Y]f 



(98) 



where Xj' is the transpose of X_. We use the fact that Xs and X^c are independent to expand 
^ as 

E^ 

E 



var [ Y,Xk\XrXs^ 

k&S 



var ( ^Xk 

ykeS 



Yl ^kXr^S 

ikes 



Pr 



(99) 



where X^.^ = (X^ — E{Xr \ X^c)) is a Gaussian random variable with variance 
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Pr,S = E 



X„ 



E[var{Xr\Xsc)]. 



(100) 



Substituting (|99l) in (|93]) and using ([5]) to bound var (X^) for all k, we obtain, 



/ 



i?.s<C 



Y^ var (Xfc) p 
fee-s 



1 Z72 



,5 



E^ 



fce5 



\ 



<c 



x^ 

V 



,fc65 



X, 



(101) 



(102) 



We define 'jk, for all /c G /C, by 



E [XfcX,] ^ x/^^Tv^. 



(103) 



Note that by definition, 



7fc G [0, 1] for aWkeJC 



(104) 



and 



K 



5^7. <1. 



fc=i 



Using the independence of X^ for all k and (|103l) . we write 



E 



2jXfcXr 



.fcG5 



E^[XfcX,] = EV7fenn. 



kes 



kes 



Next we use (11031) to evaluate Pr,5- We start by considering the random variable 



(105) 



(106) 



Xr — X^. — E [X.r|X/fJ . 



(107) 



Using (11031 ) and the independence of X^ for all k, we can write the variance of X^ as 



E 



X„ 



E[mr(X^|Xi^) 

(1 - 7^) Pr. 



(108) 
(109) 
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where we used (|98]) to simplify (11081) to (I109I ). Continuing thus, we consider the random variable 



v\fi — Ji^r — h/ 



Xr\X 



r\^K-l 



. Using the independence of Xk for all k, we thus have 



E [X^] = E 



X^ 



E 



E' 



XrlXx-l 



E[var{Xr\XK-iXK)] 

(1 -7a'-1 -lK)Pr- 



(110) 

(111) 
(112) 



Generalizing the above, we have 



1 - Y.lk]Pr = IscPr for all S C)C. 

keS" 



E [var {Xr\Xsc)] -- 
Finally, we substitute (11131 ) and (11061 ) in (IIOII ) to simplify the first bound as 



(113) 



kes 



c{E^.]. ifE7. = l 



Rs< < 



/ 



c 



V 



E 

k€S 



Pk, 

Nr 



fces 



'\ 



/ 



fcG5<= 



Otherwise. 



(114) 



Observe that for K = 1, we have Vi = Xr and 71 = 1, and thus, (ITOl) simplifies to the first 
outer bound in [3, theorem 5] for the classic single source degraded relay channel. Finally, from 
(11131) . observe that 7^., for all k, satisfies 



Eik<i. 

kf^K. 

Consider the bound B^s in ^ with U a constant. Expanding Y^ using ^, we have 



(115) 



Rs<C E 



var ( '^Xk + Xr\Xs^ 
kes 



N, 



C 



Z[Pk + 2E{ XkXr,s ]]+E [var{Xr\Xsc 

kes 



N, 



V 



J 



(116) 



(117) 



Using (El), (fm ) and (fT06l) . we simplify (fTTTl) as 



Rs<C 



Y.Pk + ls^Pr + '^T.VW^r 

ke_s kes 



(118) 
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Writing 5^,5 and B^^s to denote the bounds on the right-side of (|114l) and (11181) . respectively, 
we have for a constant U, 

Rs < min {Br,s, Bd,s) for all 5 C /C. (1 19) 

Appendix II 
Inner and Outer Bounds: Polymatroids 

We first prove that the rate regions TZ°^ and 71°]^ given by the cutset bounds are polymatroids. 
Using similar techniques, we then show that the DF regions IZ,. and IZd are polymatroids. 

A. Outer Bounds 

Consider the set functions (see [ST]) 

[ 5 = 

and 

f / (Xs; YJXscXrU) 5 C /C, 5 ^ 

[ 5 = 

for some distribution satisfying ([8]). We claim that /i and /2 are submodular [15, Ch. 44]. To 
see this, we first consider /i and ki, k2 in /C with ki ^ k2, ki ^ S, k2 ^ S, and expand 

MS U {k,}) + /i(5 U {k2}) = I{XsXk,Xr; F,|X(5u{fc,})cf/) + /(X^X^.X,; F,|X(5u{fc,})cf/) 

(122) 

= J(X,,; F,|X(5u{fci})cf/) + I{XsXr; Yd\XscU) (123) 

+ IiXsXk,Xr; F,|X(5u(fc,})ct/) (124) 

where (|123l) follows from the chain rule for mutual information. We lower bound the first term 
in (fT23]) as 

/i(XfcjX(5u{;,j)cf/) - /i(XfcjX(5u{fcj)cFdt/) (125) 

= /i(XfcjX(5Lj{fc^^fc2})cf/) - /i(XfcjX(5u{fci})cFdf/) (126) 

> I{Xki;Yd\X(^su{ki,k2})'^^) (12'7) 
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where (11251) follows from the Markov chain Xj, - U - Xj for all kJelC,k^ j and (11271) 
because conditioning cannot increase entropy. The expression (11271) added to the final term in 
(fT23])is 

Inserting (|127l) into (|123l) . we have 

MS u {h}) + /i(5 u {^2}) > MS) + MSU {fci, h}) 

for all 5 C /C. The set function /i(-) is therefore submodular by [15, Theorem 44.1, p. 767]. 
The above steps show that the rate region 7^^^ defined by the destination cutset bounds (see 

©) 

Rs < HXsXr; Yd\XsM), SCJC (128) 

is a polymatroid associated with /i(-) (see [15, p. 767]). 

One can similarly show that /2(-) is submodular. To see this, consider /2 and ki, ^2 in K, with 

ki ^ k2, ki ^ S, k2 ^ S, and expand 

f2isu{k,}) + Msu{k2}) 

= I{XsXk^;Yr\X(^s^j!^i,^y^cXrU) + li^XsXk^] Yr\X(^svj{k2})C XrU) (129) 

= I{Xkx]Yr\X(^syj{ki})cXrU) + I{Xs] Yr\XscXrU) + /(X^Xfcj ; y;|X(5u{fc2})cXrf/) (130) 

where (|130l) follows from the chain rule for mutual information. We lower bound the first term 
in ([T301) as 

h{Xk^\X^^SyJi^^.^^^)C XrU) — h{Xk^\X(^sii{k^})cYrXrU) 

= /l(XfcJX(5u|^^^fc2})'^-^rf^) — ^(-^fci|-^(5U{A:i})C^r-^r^) (131) 

> IiXk-t;Yr\X(^S^{ki,k2})'^^rU) (132) 

where (11311) follows from the independence of Xk and (11321) because conditioning cannot increase 
entropy. The expression (11321) added to the final term in (11301) is 

IiXsu{kl,k2}'^r\X(^S^J{kl,k2})(^XrU). 



35 

Inserting (|127l) into (11231) . we have 

f2iS U {h}) + /2(5 U {^2}) > f2iS) + /2(5 U {h, k2}) 

for all S C }C. The set function /2(-) is therefore submodular by [15, Theorem 44.1, p. 767]. 
This in turn implies that the rate region 7^°* defined by the relay cutset bounds (see (17])) 

Rs < I{Xs; Yr\Xsr^XrU), SCK. (133) 

is a polymatroid associated with /2(-) (see [15, p. 767]). 

B. Inner Bounds 

For the inner DF bounds, we consider the set functions (see [5T]) 

[ 5 = 

and 

, .^,_\ nXs;Yr\XscV^XrU) 5C/C,5^0 

[ 5 = 

for some distribution satisfying (|52l) . The functions /i(-) and /2(-) differ from /3(-) and fi{-), 
respectively, in the absence of the auxiliary random variables Vjc- The proof of sub-modularity 
of /a and f^ follows along the same lines as those for the outer bounds except now we have the 
Markov chain (X^, 14) - f/ - {Xj, Vj) for all k ^ j. 

We thus have that the rate region TZr defined by the DF relay bounds (see (|5T]) ) 

Rs < I{XsXr; YdlXscVsc), S C IC (136) 

is a polymatroid associated with /3(-) (see [15, p. 767]). Similarly, the region 7^^ defined by the 
DF destination bounds (see (|5TI) ) 



Rs < I{Xs; Yr\Xs^V,cXr), 5 C /C (137) 
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is a polymatroid associated with f^^-) (see [15, p. 767]). 

Appendix III 
Concavity of Bd,s and /^ .5 
A. Outer Bound B^s 

Recall that the cutset bound at the destination, Bd,s, is given by 



Y.Pk I 1- E -Tfc I -P- 2 X; VlkPkPr 



Bd,s = C 



kes 



keS" 



N, 



+ 



fce5 



N. 



V 



for all S CIC. 



(138) 



/ 



We show that Bd,s is a concave function of 7 . To prove concavity, one has to show that the 



Hessian or second derivative of -8^,5, V B^^s, is negative semi-definite, i.e, x V Bd,sx < for 



all ^ G 7^^ [21, 3.1.4]. We write 



Bd,s = UoglKo + 2Y, KkVl^] 

V k€S J 



(139) 



where 



E Pk 



TV- _ 1 I kes I Pr(l-c) 
Ao - i + ^7— + ^v^ 



T<r. — / Pk P'- 



keS. 



(140) 



The gradient VBd^s is given by 



VBd,s = [dBd,s/d^k]keK 

1 r iT 



1 



Vs Us- 



(141) 
(142) 

(143) 



where v^ is an |5| -length vector with entries Vk = K^ /^/Tfe for all k E S, v^c is an |5'^| -length 
vector with entries v^ = — P^ /Nd for all m E S^, and 



\ kes J 



(144) 
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The Hessian of 5^,5, V'^B^^s, is given by 



V'5,,5 = [d'Bd,s/djkdjm] 



\/k,m(^K 



where 



^diagi 



72 (V5, 



zz 



(145) 
(146) 



d,S) 



ds dg'^ 



(147) 
(148) 



such that d^ is an |5| -length vector with entries d^ = K^ / '^IfJ for all k E S, and dgc is an 
I iS'' I -length vector with entries dk = —2P^ /{NJKg) for all k G S"". Using the fact that Kk and 
7jk are non-negative for all k, from (11461) . for any x E TZ^ , we have 






<0 



(149) 
(150) 



with equality if and only if a; = 0. In proving the concavity of -8^,5, we assume only that 7^ > 0, 

for all k. Thus, from continuity, the concavity also holds for all non-negative 7^ satisfying (see 

dn])) 

Y.^k<l. (151) 

fce/c 



For a fixed 7 , we now find the 7 that maximizes Bds subject to (|151l) above. For a 

——C> ——C> ' 

c G [0, 1), we fix 7 such that its entries 7^, for all k E S'^, satisfy 



and thus, from (|151|) we have 



^ 7fc = 1 - c, 



X^ 7fc < c. 
fce5 



(152) 



(153) 



Since Bd,s is a continuous concave function of 7 it achieves its maximum at a 7 J where 






for all A; G 5. 



Tfc 



(154) 
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Using the method of Lagrange mukipliers, we find that a 7J that maximizes i?^^ subject to 
(11521) and (11531) has entries 

7^ = I T% ^^^ ■ (155) 

l^ kes 

B. Inner Bound Id,s 

Recall that the DF bound, Id,s, at the destination is given as 

Y.Pk U-Y.Pk]Pr 2 X] ^/{l~Cik)l3kPkPr 



i-d,S 



c 






+ 



fc65= 



kes 



Nr, 



K, 



V 



for all S CIC. 



(156) 



/ 



Comparing (11381) and (11561) . for 7^ = (1 — a^) P^ for all A; G 5 and 7^ = (3^ for all k G S'^, the 
DF rate bound in (11561) simplifies to that for the outer bound in (11381) . and thus, one can use the 
same technique to show that 1^,5 is a concave function of a^^ and P . For the power fractions 
/3k, we have 

5]/?fe<l- (15V) 



keK 



For a fixed a^^, we determine the optimal P maximizing 1^,^ by fixing the vector /5 such that 

J2p, = 1-c (158) 



fe65= 



J]/3fc<c. 



(159) 



kes 



where c e [0, 1). Since Id,s is independent of /S^ for as = 1, we assume that as 7^ 1- 

We now consider the special case in which a^ ^ I and /3 are fixed. We determine a /5 
that maximizes Id^s subject to (11591) and (11581) . Since Id,s is a continuous concave function of 
/? it achieves its maximum at a /JJ where 



dld,s 



aPk 



Pt 



for all keS. 



(160) 



As before, using Lagrange multipliers, the optimal /?J that maximizes 7^5, subject to (11591) . has 
entries 

« = |if7Sa^ ^^•S . (161) 

Rate region for a fixed ox- For any choice of a non-zero Ojc and a Pjc satisfying (11571 ). the 
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(J 



Pi 
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P^/N^ = P^/N^ = 7 dB 
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la Id r d 
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- -Q- (.15,.85) 
■■■+■■■ (0,1 
^^^— Sum Bound 







Rate R (bits/ch. use) 




Fig. 3. Rate region achieved at tlie destination for a two-user MARC and Ofi = a2 = 1/2. 



rate region given by (11561) for all 5 is a polymatroid. For ojc = \., from (11561) we see that there 
are no gains achieved from coherent combining, i.e., it suffices to choose Ptc = 0. Consider 
Ox: 7^ 1. Since there is at least one k for which a^ < 1, gains from coherent combining at the 
destination are maximized by choosing (3}c to satisfy (11571) with equality. For a fixed ox, we 
then write the rate region at the destination as a union over all polymatroids, one for each choice 
of (3]c satisfying 

~ K 

5^/3, = 1. (162) 



fc=i 



Observe that for 13* with entries given by (11611) . the bound Id,s is maximized. In Fig [3l we 
illustrate the rate region for a two-user degraded Gaussian MARC with the SNR Pi/Nd = P2/Nd 
chosen as —10 dB, a = (1/2, 1/2), and five choices of Pic- Observe that the maximum single- 
user rate -Ri is achieved by setting Pi to 1 though this value does not maximize R2 or Ri + R2. 
For all other (/3i,/52) such as (0.85,0.15), as /3i decreases and P2 increases, Ri decreases while 
i?2 increases achieving its maximum at /32 = 1- The bound on the sum rate Ri + R2 increases 
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from (/3i,/52) = (1,0), achieves its maximum at {Pl,P2) — (V^, 1/2), and then decreases as 
P2 approaches 1. The resulting region at the destination is shown in Fig. [3] as a union over all 
polymatroids, one for each choice of Pic. 



Appendix IV 
Br,s vs. J_ic 
We show that the function B^ s in (flOl) is a concave function of 7 for a fixed 7 and for 
all iS C /C. Recall the expression for 5^,5 as 



/ 



Br,S = C 



Pk 



E VikPk 

kes 



\ 



E 

\ \ fce5^ / / 



(163) 



where we assume that 



^ 7fc = 1 - c < 1. 



(164) 



keS'^ 



Observe that 5, .5 is maximized when c = 1, i.e., 7a; = for all k G S, and minimized for c = 0. 
Further, comparing Br^s and Bd,s, one can see that for 



7fc 



Pk/{EkesPk) , keS 
, keS^ 



(165) 



Bf^s achieves its minimum, i.e., B^^s = 0. 
We write 

^VikXk 



X 



Kkes 



where 



Pmax = maxfcex; Pk and A^ = Pk/ Pn 



(166) 



(167) 



Substituting (|166l) in the expression for 5^,5 in (|163l) . we have 



Br 



^ I sr~^ k ■^ -'max 



(168) 
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Differentiating B^^s with respect to x we have 



dBr,s 


—rP ( Pi r^P 

■^-' max ^ 1 Y^ '^ max 


dx 


/ \ 2 


d Br^S 


N.C y' 1 ^^^^. j (, ^.e ; 


dx^ 


/ \2 



(169) 



(170) 

\ k&S J 

<0 (171) 

where the strict inequality in (11711) follows since all terms in (11701) are positive. Further, for any 
c > , from (|168l) i?,. 5 is maximized at x = 0, i.e., for 7a; = for all k E S. Thus, we see that 
Br^s is a concave decreasing function of x. 

Appendix V 
Proof of Theorem [5] 

We now prove Theorem |5] and give the solution to the max-min optimization 

^^ = ."If^, "^^^ {^^'^ (lye) ' ^^'^ (lyc) ) • (1^2) 



Consider the function 



J(7^, 6) = 6Br,,c (7^) + (1 - ^) 5d,yc (7^) , 6e[0,l]. (173) 

Observe that J is linear in 6 ranging in value from J^x: for (5 = to Ir^jc for 6=1. Thus, the 
optimization in (11731) is equivalent to maximizing the minimum of the two end points of the 
line J over Tqb- Maximizing J(7„,(^) over 7 , we obtain a continuous convex function 

V{5)= max J{^^,5), <5e[0,l]. (174) 

From (11731 ) and (11741 ), we see that for any 7 , ^(7„, 5) either lies strictly below or is tangential 
to V{5). The following proposition summarizes a well-known solution to the max-min problem 
in (flTll) (see [9]). 
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Br.M 




B,x(rK) 



KM 




KM 



hM 




5 



hM 



Fig. 4. Illustration of Cases 1, 2, and 3. 



Proposition 3: 7* is a max-min rule where 



5* = arg min V(6). 

<5e[o,i] 



(175) 



The maximum bound on Rjc, V{5*), is completely determined by the following three cases (see 
Fig. HI). 



Case 1: 6* = : ¥(5*) = 5,,^(7* J < BrAl*) 



'-JCS* 



Case 2: 5* = 1 : ¥(6*) = 5.,^(7^,J < B.^^ij*^ J 



LJCS* 



Case 3: < 6* < 1 : ¥(6*) = 5,,yc(7^ ,J = Bd,icir^ 



K,&* 



(176) 

(177) 
(178) 



We apply Proposition [3] to determine the maximum bound on /?/c. We study each case 
separately and determine the max-min rule 7* for each case. In general, the max-min rule 7* 
depends on an optimal 5* . However, for notational convenience we henceforth omit the subscript 
5* in denoting the max-min rule. We develop the optimal 7* and the maximum sum-rate for 
each case. We first consider case 1 and show that this case is not feasible. 

Case 1: This case occurs when the maximum bound achievable at the destination is smaller 
than the bound at the relay. In Appendix Hill we show that the bound Bd,K.{l^) is ^ concave 
function of 7 and achieves a maximum at Yk, whose entries 7^ satisfy (I115h and are given as 



'-K 



ll = Pk/{Y.keicPk), forall/cG/C. 



(179) 



Substituting (|179l) in (flOl) . we have Br,K{lt) = which contradicts the assumption in (11761) . 



^JJC 
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thus making this case infeasible. 

Case 2: Consider the condition for case 2 in (11771) . This condition implies that the case occurs 
when the maximum bound achievable at the relay is smaller than the bound at the destination. 
From (l59l) . we observe that Bj-f^ decreases monotonically with 7^ for all k and achieves a 
maximum of 



i?.,^(7y=C[^ (180) 

^t 7* = 0. Comparing (flOl ) and (fTTl) at 7* = 0, we obtain the condition for this case as 

Case 3: Finally, consider the condition for Case 3 in (|178l) . This case occurs when the maximum 
rate bound achievable at the relay and destination are equal. The max-min solution for this case 
is obtained by considering two sub-cases. The first is the relatively straightforward sub-case 
where 7* = is the max-min rule. The resulting maximum sum-rate is the same as that for 
case 2 with the condition in (|181l) satisfied with equality. Consider the second sub-case where 
7*. 7^ 0, i.e., when 

Yl,Pk Yl,Pk + Pr 

We formulate the optimization problem for this case as 

maximize Br^K (7) 

subject to B,.ic (7) = Bd,ic (7) ■ 

We write 

-Pmax = maXfcgx: -Pfc) ^k = Pk/ Pmax, (184) 

and define 



(183) 



^=./EMfc. (185) 

fce/c 
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Substituting (11841) and (11851) in ([JO]) and ([IT]), we have 



ff 



Br,d^) = C 



Vfcg/c / 



(186) 



/ 



5d,yc(x)=C 



kaK. 



N. 



\ 



(187) 



/ 



Observe that Br^/cix) and Bd,ic{x) are monotonically decreasing and increasing functions of x, 
respectively, and thus, the maximization in (11831) simplifies to determining an x such that 



7 J -'k -^ -'max 7 y -'A; T -f^r ~r -^S^V -'max-')" 

fce/C fc6/C 



We write 



Nr 



Kq — -Pmax /Nr 



Nd 



i^2 = %^ + #-, and i^3 



Ki = vK:j^-/Nd 



(188) 



(189) 



Nd ' Afd' """ ""■^ Nr 

From (f82l) . since i^3 > 7^2, the quadratic equation in (|188l) has only one positive solution given 

by 



X 



K, 







The optimal power policy for this case is then the set Q of 7* for which 7* satisfies (11851) with 



-Lie 



-Lie 



X = x* in (11901) . The maximum achievable sum-rate for this case is then obtained from (11861) as 

(191) 



C 



keic 



Nr 



J 



Appendix VI 
Proof of Theorem [TTl 

We now prove Theorem \TT\ and give the solution to the max-min optimization 



R,c = max min ( I,,^ (a,^) , J^,^ ( a^, p )) 



{^K^^^)'^^ 



(192) 



As in Appendix |Vl a solution to the max-min optimization in (11921) simplifies to three mutually 
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exclusive cases [16, II.C] such that the max-min rule (aj^,/?* ) satisfies the conditions for one 
of three cases. The conditions for the three cases are 

Casel: Id,jc{(^h Sc) < ^r,jc{(^l) (193) 

Case 2: Ir,jc{(^l) < Id,jc{ah (T^) (194) 

Cases: /.,yc(«^) = /d,yc(«^,^;^) (195) 

We develop the conditions and determine the max-min rule for each case. We first consider case 
1 and show that this case is not feasible. 

Case 1: This case occurs when the maximum bound achievable at the destination is smaller 
than the bound at the relay. Observe that Id,K{(^K^ Pir) ™ <l60l) decreases monotonically with a^, 
for all k, and, for any [3tc, achieves a maximum at aj^ = of 

Id,d<y.hP^) = C — j;^ I . (196) 



However, substituting a*^ = Oin (1591 ). we obtain 

/r,yc(«^) = (197) 

which contradicts the assumption in (11931 ), thus making this case infeasible. 

Case 2: Consider the condition for Case 2 in (11941) . This condition implies that the case 
occurs when the maximum bound achievable at the relay is smaller than the bound at the 
destination. From (|59| ), we observe that J^^/c increases monotonically with at for all k and 
achieves a maximum of ^_ 

/EPk\ 

IrA^l) = C ^ (198) 

at ajc — i- Comparing (|59| ) and (l60l) at aj^ = \, we obtain the condition for this case as 



X] -Pfc YliPk + Pr 

^^'^ < ^^^ (199) 

Nr - Nd - 

Case 3: Finally, consider Case 3 in (11951 ). This case occurs when the maximum rate bound 
achievable at the relay and destination are equal. The max-min solution for this case is obtained 
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by considering two sub-cases. The first is the relatively straightforward sub-case where ajc ~ 1 
is the max-min rule. The resulting maximum sum-rate is the same as that for case 2 with the 
condition in (11991) satisfied with equality. Consider the second sub-case where aj^ ^ \, i.e., 

k^j^_ k&K ^200) 

In Appendix Hill we show that, for a fixed Ojc, Id,s, is a concave function of f3jc for all S C JC. 
Furthermore, from (l57l) . for a;c 7^ i> Id,K in (|60l) is maximized by a /3* whose entries Pi, for 
all k e IC, satisfy 

EPl = ^ (201) 

fceK: 

and are given as 

p* = <j ELi(i-"fe)Pfc -/c ^ - ^^^ ^^^ ^ ^ ^^ ^202) 

Observe that the optimal power fraction Pi that the relay allocates to cooperating with user k 
is proportional to the power allocated by user k to achieve coherent combining gains at the 
destination. Thus, one can formulate the optimization problem for this case as 

maximize Ir,jc («) 

subject to Irx (a) = h,K. («, ^) , (203) 

Using Lagrange multipliers we can show that it suffices to consider P^ = Pi in the maximization. 
Since the optimal PI in (|202|) is a function of ox, Id,K:{(^K:->I^V) simplifies to a function of a^^ 
as 



h,ic{^K,P*)=C 



K' 



(Y.Pk + Pr + 2 Y. il-ak)PkPr\ 
keK. \ kGK. 



N. 



V 



(204) 



We further simplify Id,K{Qi_fQ,P*) and Irxif^jc) ^^ follows. We write 

-Pmax = maXfcg/c Pk, ^k = Pk/ Pma.x, (205) 
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and 



1= JH (l-«fc)Afc. 

keK 



(206) 



Substituting (12051) and (12061) in ^9^ and (160]), we have 



// 



Ir,M = C 



\keK J 






\ 



Nr 



(207) 



ff 



W{q) = C 



(EPk]+Pr + 2qVK:^r 
\keK J 



\ 



Nh. 



V 



(208) 



/ 



Observe that Ir,ic{Q) and Id,K{<l) ^^ monotonically increasing and decreasing functions of q and 
thus, the maximization in (12031) simplifies to determining a q such that 



EPk- q'Pr 
k&K. 



max / J -^ fc ~r -fr ~r ■^Q'V -'max-'r 

fce/c 



(209) 



The condition in (12091 ) has the geometric interpretation that the bounds on i?^: ^e maximized 
when the Ji'-user sum rate plane achieved at the relay is tangential to the concave sum-rate 
surface achieved at the destination at its maximum value. We further simplify (12091) by using the 
definitions in Appendix IVl for Kq, Ki, K2, and K3. From (12001) . since K3 > K2, the quadratic 
equation in (12091) has only one positive solution given by 






(210) 



The optimal power policy for this case is then the set V of {a^, P*{a'^)) such that a*^ satisfies 
(12061 ) for q = q* and for each such choice of a*^, (3* is given by (12021 ). The maximum achievable 
sum-rate for this case is then given by 



C 



^ ^Y.P^-{q*?Pnr.} 
keK. J 



Nr 



\ 



(211) 



/ 



Remark 12: The optimal q* in (12101) is the same as the optimal x* in (11901 ). Further, the 
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maximum inner (DF) and outer bounds on the sum-rate are also the same for the equal-bounds 
case in (12111) and (11911 ). respectively. 



Appendix VII 
Sum-Capacity Proof for the Active Class 

In Theorem [121 we proved that DF achieves the sum-capacity for an active class of degraded 
Gaussian MARCs. In the proof we argue that since the maximum DF sum-rate is the same as 
the maximum outer bound sum-rate, every DF max-min rule (aJc'/^t-) ^ ^a '^hat achieves this 
maximum sum-rate, i.e., for which Rria}^) H Rd{(^K.^ Pt) belongs to the set of active cases, also 
achieves the sum-capacity. We now present a more detailed proof of the argument. 

We begin by comparing the inner and outer bounds. As in the symmetric case, without loss 
of generality, we write 

7fc = (l-afc)/3fc for all A; (212) 



where {aj^, (3 ) G T. We then have. 



"-IC 



( 



5r,^(«yC,/5J=C 



K.' 



Pi 



E 



k \k€JC 



Oik) PkPk 



'■\ 



Nr 



(213) 



/ 



and 



J2Pk + Pr + 2j:y/{l-ak)PkPkPr\ 



B,A^!C,Pr)=C 



fce/c 



fce/c 



■K^ 



N. 



h,K.{.Q^K.i (^ 



K> 



(214) 



/ 



Choosing /5 as the DF max-min rule /?* in (TtTI) . simplifies (12131) to 



BrA^K.i^l)=C[Y. 



-K. 



C^kPk 



ItA^k)- 



(215) 



Using theorem [TTl one can then verify that -Br,/c(a5c' P*jr) = Bd,K:{(y*K:-, (^V) i^ achieved by all o*^ E 
V. Consider a. a^ E Va and a corresponding (3^ such that the DF region T^rlajc) ^ T^di^lK., PV) 
belongs to the set of active cases. From Theorem [TTl this implies that 



Id,Ai^hP*r) + IrM^h) >P = B* for all AC)C. 



-K 



(216) 
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Using (|212l) . we expand i?^^ in (fTTI) as a function of (aj^,/?* ) as 

/ 



B,^s{ak,P*)=C 



K' 






A^^ 



/ 



(217) 

>/d,5(a^,^;) (218) 

where (12181) follows from the fact that (1 - al) (31 < (31, for all k and for all (a^,^)^). It is, 
however, not easy to compare -Br,5(ttA;'/^I-) ^^^^ Ir,s{(^K)- Note, however, that the choice of 7^ 
in (|212l) requires the same source-relay correlation values for both the inner and outer bounds. 
Furthermore, for every choice of Gaussian input distribution with the same K correlation values 
for both bounds, comparing the degraded cutset and DF bounds in ^ and (|5T1) . respectively, 
for a constant U, we have 

I{Xs; Yr\XscXr) > I{Xs; F.lX^cl/^X,) for all 5 C /C (219) 

where in (12191) we use the fact that conditioning does not increase entropy to show that the 
cutset bounds at the relay are less restrictive than the corresponding DF bounds. From (|215l) . 
the inequality in (12191) simplifies to an equality for S = JC and for (a*^, P*) E Va when 7^ is 
given by (12121) . Combining these inequalities with (|216l) . we then have 

BdAi^lhSc^ + BrMo*^) >I* = B* for all A C JC. (220) 

Thus, every DF max-min rule that results in an active case polymatroid intersection, i.e., every 
(f^Jc'/^ir) ^ ^d' ^1^0 results in an active case for the outer bounds when 7^ is given by (12121) . 
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